Strong CP breaking and quark-antiquark repulsion in QCD, at finite 9 
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This work is devoted to the study of the CP- breaking dynamics in QCD, at finite S-angle. By 
working in the semi-classical limit, in which the topology of the vacuum is clustered around in- 
stantons and anti-instantons, we show that the quantum fluctuations of the ^-vacuum generate 
an effective flavor-dependent repulsion between matter and anti-matter. As a consequence, during 
the tunneling between the degenerate classical vacua, quarks and anti-quarks in a neutron migrate 
in opposite directions, giving rise to an oscillating electric dipole moment. We discuss a possible 
phenomenological implication of this effect. 



The hypothesis of CP-violation in the strong sector 
of the Standard Model arises from the analysis of the 
quantum structure of the QCD vacuum. At a classical 
level, the QCD Hamiltonian possesses an infinite set of 
degenerate gauge-dependent vacua, each of which is char- 
acterized by its winding number n. The gauge-invariant 
definition of the ground-state is given in terms of a linear 
superposition of such classical vacua (9- vacuum) PJ: 
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The structure of the QCD vacuum, combined with the 
observation of weak CP-violation in neutral kaon sys- 
tems, gives rise to the so-called ^-term in the QCD action 
(in the Euclidean formulation): 
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where 8 = 9 + argdet(M), and M is the complex, non- 
hermitian quark mass matrix, arising from the sponta- 
neous breaking of the electro- weak gauge symmetry. The 
real constant 9 is an additional parameter of the Stan- 
dard Model, which has to be fixed by experiment. At the 
moment, the most constraining bounds on this quantity 
come from measurements of the neutron electric dipole 
moment (EDM) @, which indicate that 6 < 1CT 9 0,0. 
In this Rapid Communication we explore some dynami- 
cal consequences of the 0-term at the microscopic level. 
In particular, we identify the mechanisms which gives rise 
to a neutron EDM in the semiclassical limit of QCD, in 
which the mixing between the degenerate vacua in Q is 
mediated by instantons. 

Our main result is the discovery that, during the tun- 
neling processes, the #-term (|2J) generates an effective 
repulsion between matter and anti-matter, in the neu- 
tron. As a consequence, quarks and anti-quarks migrate 
in opposite directions giving rise to a finite EDM. Our 
results show that, at least on the semi-classical level, the 
EDM arises from the local separation of positive and neg- 
ative baryonic charges in the neutron. It docs not follow 
from the displacement of the positive and negative elec- 
tric charge carried by the valence quarks, as one would 



intuitively expect in a naive non-relativistic quark model 
picture. These results are obtained comparing the den- 
sity distributions of the baryonic and electric charges. 
Our starting point is the neutron baryon number den- 
sity correlator, defined as Gb(y, t) = 1/3 Gd(y, t) + 
1/3 G„(y, r), where 
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In © t is the Euclidean time, Jt{x) = Qf, ^3 is a 
diagonal matrix of spinor indices T3 = diag(l, — 1, 1, —1) 
and Jn(x) = £abc(wa G75 db) d c is an interpolating field, 
which excites states with the quantum numbers of the 
neutron. Similarly, we define the electric charge density 
correlator as G e/m (y,r) = 2/3 G„(y,r) - 1/3 G d (y,r). 
The correlation function @ measures the probability 
amplitude to find a quark of flavor / at the point y, 
in a system with neutron quantum numbers. 

The component of the EDM along the direction of the 
neutron spin can be extracted from the Green's function: 
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which, in the large time r limit, reads: 
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where D z is the component of the neutron EDM along 
the spin direction, M is the neutron mass and A^r is the 
coupling of the neutron to the interpolating operator, 
(0\J N (0)\N) = A N u p . 

In QCD, a finite dipole can only arise from the CP- 
breaking interaction |2|l. To lowest order in 8, we can 
write (in an obvious notation): 



G/(y,- 
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In the semi-classical limit, the topological charge is 
condensed around instantons and anti-instantons (Q = 
Ni — Na), and the physics of the quantum mixing of 



2 



Q=-28 



at 

IE 



y 3 [ fm l 
(A) 



10 



G of = d 



Q=0 



TiTTTT 

i it 



y 3 [ fm l 
(B) 



Q=+28 



CD 



■ ■ f = u 

o of = d 



y 3 [fm] 



(C) 



FIG. 1: The baryon number density correlators for u and d quarks, Eq. in units of 6 x 10 5 fm 12 , calculated in canonical 
ensembles with topological charge Q — —28, Q — 0, Q — +28. 



the 0-vacuum can be formulated in terms of an intuitive 
pseudo-particle picture. The path-integral can then be 
computed by summing over the configurations of a sta- 
tistical grand-canonical ensemble of pseudo-particles. In 
such an approach, a quantitative estimate of the matrix 
element (JHJ can only be performed in a model-dependent 
way, as we do not know from first principles the den- 
sity and size distribution of pseudo-particles in the vac- 
uum. This is the starting point of the Instanton Liquid 
Model (ILM) - for a review see 8j - which has been proved 
to be very successful in describing the phenomenology of 
the QCD vacuum and of light hadrons. 

In the present study, we choose to avoid introducing 
model-dependent parameters and we focus on the qual- 
itative effects generated by the dynamical interplay of 
the strong CP-breaking interaction (j2J with the quan- 
tum structure of the 9- vacuum, at semi-classical level. A 
quantitative, albeit model-dependent, prediction of the 
neutron EDM in this model will be presented in a sepa- 
rate publication . 

In the instanton vacuum, matrix elements with one in- 
sertion of the topological charge operator can be written 
as a sum over the contributions of the different topolog- 
ical sectors [Toj | : 

(°3^a/ ^W)=EWI)<J(0> 9 , (7) 

Q 

where V(Q) denotes the relative occurrence of configu- 
rations with topological charge Q and (0)q is the aver- 
age performed in a canonical ensemble with such a total 
topological charge. This equation expresses the fact that 
CP-breaking interactions arise from topologically non- 
trivial gauge configurations. Physically, it implies that 
these forces are triggered only during tunneling between 
distinct classical vacua. 

As long as we are interested in qualitative phenom- 
ena, we do not need detailed knowledge of the positive- 
definite weight factor PflQI) in Q. We shall there- 



fore concentrate on the contribution from each topolog- 
ical sector (factor Q (0)q in Q). This can be done by 
evaluating averages in different canonical ensembles in 
which the number of instantons and anti-instantons is 
fixed. In order to do so, we have used the Interact- 
ing Instanton Liquid Model (IILM), developed in [llj . 
We have averaged over 1000 configurations of an ensem- 
ble of 130 pseudo-particles in a periodic box of volume 
V = 3.2 x 4 fm 4 . For example, the contribution of 
the topological sector Q = 10 was obtained inserting 70 
intantons and 60 anti-instantons in the box. For each 
configuration, the fermionic determinant and the quark 
propagator have been calculated by diagonalizing numer- 
ically the Dirac operator (for a detailed discussion of 
the method see Section VI in 8]). In topologically non- 
trivial sectors, the quark propagator receives contribu- 
tion from -/Vf |Q| exact zero- modes (semi-classical realiza- 
tion of the Index Theorem) . The ensemble of collective 
coordinates of the pseudo-particles was then generated 
dynamically, through an usual accept/reject Metropolis 
algorithm. To improve the signal-to-noise ratio, we have 
used rather large quark masses, m u = = 75 MeV 
and m s = 150 MeV, and a rather small Euclidean time, 
t = 0.7 fm. We note that, in the present approach, we 
do not make use of an anomalous Ward identity to re- 
late the topological charge to the divergence of a gauge- 
invariantly defined axial current. Hence, the results of 
the present EDM calculation do not allow to distinguish a 
spontaneous from an anomalous [j| instanton-induced 
breaking of the axial symmetry (for an example in which 
the two types of breaking lead to different predictions, 
see @). 

In Fig. we report the results of the averages of the 
baryon density correlators © performed in canonical 
ensembles with topological charge Q — —28, Q = 0, 
Q = +28 ( with y chosen along the z direction ). Cal- 
culations in ensembles with other positive and negative 
total topological charge have also been performed and 
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FIG. 2: Contribution of the topological sectors Q — —28 and Q = +28 to the neutron baryon number density correlator (left 
panel) and to the neutron electric charge density correlator (right panel), in units of 8 x 1CP 5 fm~ 12 . A negative baryon 
number density for j/3 > denotes a local accumulation of anti-matter in the positive yz hemisphere. 



the results were found to follow the same trend. Hence, 
at a qualitative level, Fig. QJV (Fig. CP) can be regarded 
as representing all terms (0}q in with Q < (Q > 0). 
(We choose to plot the results of the topological sectors 
with a large topological charge, because they have a bet- 
ter signal-to-noise ratio.) 

These results have several interesting implications. 
First of all, we note that the baryon density at zero net 
topology ( Fig. 03 ) is even under parity transformation 
y3 — * —1/3, as expected. On the other hand, contributions 
from the topologically non-trivial sectors are P-odd and 
lead to a non-vanishing EDM. 

The most interesting feature of these results is the sep- 
aration of positive and negative baryonic charges in phys- 
ical systems with neutron quantum numbers. Fig. ^\ 
and Fig. \\]P imply that in sectors with non-vanishing 
topological charge, the baryonic charge density of u 
quarks is positive for 2/3 > and negative for y% < 0. 
Conversely, the baryonic charge density of d quarks is 
negative for 2/3 > and positive for 7/3 < 0. Notice that 
this qualitative effect is the same in all positive and neg- 
ative topological sectors, due to the sign carried by the 
topological charge factor Q in (7J. 

The physical interpretation of these results is the fol- 
lowing. Any time the quantum QCD ground-state re- 
arranges itself by tunneling, the 0-term (J2J generates an 
interaction which effectively shifts the u quarks and the 
d anti-quarks along the positive direction of the neutron 
spin. At the same time, u anti-quarks and d quarks are 
shifted toward the opposite direction. The net effect is 
the creation of an electric current inside the neutron, 
pointing in the direction of the spin. 

From these correlators we can construct the contri- 
bution to the baryonic number density ( Fig. [2K )• We 
conclude that the tunneling produces a local separation 
of the baryonic charge, with an accumulation of mat- 
ter in the 2/3 < hemisphere and of anti-matter in the 
2/3 > hemisphere. Notice that the contribution to the 
total baryon charge of the neutron coming from configu- 



rations with Q 7^ vanishes, due to the odd symmetry 
of the correlators © under y% — > — j/3 transformations. 
In other words, the neutron total baryon number comes 
entirely from the topological sector with Q — 0. Hence, 
the baryon and electric charge asymmetries are associ- 
ated with the sea quarks. 

In Fig. [3)3 we show that the separation of the bary- 
onic charge induced by the 0-term generates a disentan- 
glement of positive and negative electric charge in the 
neutron. This is the microscopic dynamical mechanism 
underlying the EDM formation in QCD, at the semi- 
classical level. 

Notice that in the present discussion we have never as- 
sumed that the quarks in the correlator J2J form a bound- 
state. Indeed, these arguments hold also for short-sized 
Euclidean correlators, in which the lowest-lying pole is 
not isolated from its excitations. This means that the 
dynamical mechanism analyzed here does not only con- 
cern neutrons, but applies to all systems with nucleon 
quantum numbers. In particular, it should be effective 
also in the its excited states and in the proximity of the 
de-confinement phase transition, where the nucleon is 
melted into its partonic components. 

We insist on the fact that the outcome of our analysis 
does not depend on the particular values of the model pa- 
rameters which define the ILM. These results rely only on 
the working assumption that the quantum mixing of the 
QCD vacuum can be described in terms of isolated tun- 
neling events. On the one hand, there is no a priori guar- 
antee that a semi-classical description may lead to quan- 
titatively realistic predictions. On the other hand, such 
an approach contains the correct ingredients to point out 
explicitly the connection between the quantum structure 
of the vacuum and the effect of CP-breaking interactions. 
We believe that it is unlikely that further quantum cor- 
rections would completely destroy the qualitative mecha- 
nism which emerges from the present semi-classical anal- 
ysis. 

An interesting question to ask is how this microscopic 
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FIG. 3: Example of a process leading to a positive neutron 
EDM, in the low-energy effective description. 



QCD dynamics is translated in the language of a low- 
energy effective theory, with hadronic degrees of freedom. 
In [7j Crewther et al. performed the calculation of the 
lowest-order contribution to the EDM in chiral perturba- 
tion theory. In their approach, the EDM arises from the 
dissociation of the electric charge of the neutron associ- 
ated with the quantum fluctuation, n — * ir~p (see Fig.|2Jl. 
In order to produce a finite EDM, the neutron must dis- 
sociate into charge components, in a way that breaks the 
spherical symmetry of the electric charge distribution. In 
particular, in order to generate a positive EDM, the vir- 
tual it~ cloud must be mostly localized in the southern 
hemisphere (z < 0), relative to the direction of the spin 
of the nucleon. The same effect is achieved at the micro- 
scopic level by the non-perturbative instanton-induced 
effective repulsion discussed above. In fact, we have seen 
that the topological fluctuations drive the u and d sea 
quarks towards the region z < 0. (see Fig. G]C). 

Let us now discuss some phenomenological implica- 
tions of the non-perturbative CP-breaking dynamics dis- 
cussed in this work. One immediate consequence is 
that 8 > would imply an asymmetry in the photon- 
production of charged pions, with an excess of tt~ pro- 
duced in the z < emisphere. This effect is of the same 
type of the asymmetries discussed in |l2( , in the context 
of relativistic heavy ion collions. 

A second phenomenological implication of the present 
semi-classical description is that the neutron EDM must 
have characteristic frequency of oscillation. We have seen 
that the dissociation of the electric charge is realized 
through periodic currents, induced by local topological 
fluctuations in the vacuum. When the Q = condition is 
restored, the electric charge distribution is relaxed to its 
symmetrical equilibrium state. Hence, we expect a nat- 
ural frequency of oscillations of the EDM, of the order of 
the inverse of the topological screening lenght. Such a fre- 
quency can be estimated in different ways. For example, 
in jl.'il] we analyzed some combination of meson point-to- 
point correlators which relate directly to the amplitude 
for chirality-flips in a quark- antiquark state induced by 
topological interactions. From a spectral analysis of such 
a correlation function it was shown that the information 
about the initial chiralities of the quark and antiquark 
is exponentially destroyed in time with a characteristic 
decay-constant given by t x ~ 1/ m v > ~ \j m ao ( ao being 



the lightest 1=1, J p = + meson ). Hence, we can ar- 
gue that the characteristic semi-classical frequency of os- 
cillation of the neutron EDM is of the order 1 /m,< GeV. 

A natural consequence is that, for 8 ^ 0, we expect a 
resonance in the neutron Compton scattering cross sec- 
tion, at center of mass energy y/s ~ M + m v > ~ 2 GeV . 
The corresponding final photon will be emitted in a state 
such that the total parity will not be conserved in the 
scattering process. The observability of such a resonance 
depends on the magnitude of 8. Since the measurements 
of the EDM indicate that this parameter is at least ex- 
tremely small (if not vanishing) , it is hard to imagine that 
any sizable effect could be seen in a foreseeable scattering 
experiment. On the other hand, it would be interesting 
to study what contraints on 8 could be set from Comp- 
ton scattering data. Such a quantitative analysis would 
require a model-dependent choice of the instanton size 
and density parameters and is outside the scope of the 
present work. 

Summarizing, we have studied the dynamical mecha- 
nism for EDM formation in QCD, in the presence of a 8- 
term. We have found that the baryon number carried by 
the sea quarks in the neutron periodically undergoes a lo- 
cal rearrangement. This fact can be interpreted as due to 
a flavor-dependent quark-antiquark repulsion, triggered 
by tunneling events in the (9-vacuum. This is the mi- 
croscopic origin of the breaking the spherical symmetry 
of the electric charge distribution, at the semi-classical 
level. The same mechanism can be effective in excitations 
of the nucleon and in the proximity of the QCD phase 
transition. We have estimated the natural frequency os- 
cillation of the EDM to be of the order of the 1 /m v ' and 
argued on possible phenomenological implications. 
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